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Abstract 

A geometry-based density functional theory is presented for mixtures of hard spheres, hard 
needles and hard platelets; both the needles and the platelets are taken to be of vanishing thickness. 
Geometrical weight functions that are characteristic for each species are given and it is shown how 
convolutions of pairs of weight functions recover each Mayer bond of the ternary mixture and 
hence ensure the correct second virial expansion of the excess free energy functional. The case of 
sphere-platelet overlap relies on the same approximation as does Rosenfeld's functional for strictly 
two-dimensional hard disks. We explicitly control contributions to the excess free energy that 
are of third order in density. Analytic expressions relevant for the application of the theory to 
states with planar translational and cylindrical rotational symmetry, e.g. to describe behavior at 
planar smooth walls, are given. For binary sphere-platelet mixtures, in the appropriate limit of 
small platelet densities, the theory differs from that used in a recent treatment [L. Harnau and 
S. Dietrich, Phys. Rev. E 71, 011504 (2004)]. As a test case of our approach we consider the 
isotropic-nematic bulk transition of pure hard platelets, which we find to be weakly first order, 
with values for the coexistence densities and the nematic order parameter that compare well with 
simulation results. 

PACS numbers: 61.20.Gy, 82.70.Dd, 61.30.Cz 
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I. INTRODUCTION 



Dispersions of non-spherical colloidal particles are model systems to study a broad range 
of phenomena in condensed matter ranging from fluid phase separation to the emergence of 
liquid crystalline ordering. Examples of such systems are clay suspensions , dispersed 
gibbsite platelets P|, mixtures of silica spheres and silica-coated boehmite rods wax 
disks [5j , or nonaqueous suspensions of laponite and montmorillonite [6] . While already pure 
systems often possess complex liquid crystal phase behavior, binary mixtures additionally 
may demix into bulk phases with different chemical compositions, e.g. like in dispersions of 
disks and spheres jsj . Such mixtures are often viewed as composed of a primary species (here 
the spheres) that interact with an effective depletion potential that is generated by the second 
component referred to as the depletion agent (here the disks). The depletion interaction is 
primarily attractive with a range of attraction similar to the size of the depletant agent and 
with a strength that is ruled by the concentration of the depletion agent. Under appropriate 
conditions the depletion interaction may be sufficiently strong in order to drive a phase 
transition that is anaologous to the gas-liquid phase separation in simple substances, and 
that manifests fluid-fluid demixing when regarding the full mixture. 

Recent theoretical work has been devoted to fluids of platelike particles near a hard wall 
0], an interaction site model for lamellar colloids wetting and capillary nematization 
of binary hard-platelet and hard-rod fluids 0], the lamellar Zwanzig model 10, 11 1, and 
colloidal hard-rod fluids near geometrically structured substrates The phase diagram 
of mixtures of hard colloidal spheres and discs was obtained within a free-volume scaled 



particle approach [3l 
of platelets was studied in detai 
the Derjaguin approximation 14]. ' 



The depletion potential between two spheres immersed in a sea 
iJ . llS^ , and found to compare well with predictions from 



Sedimentation was found to influence liquid crystal phase 
transitions of colloidal platelets jl6|, as well as multi-phase equilibria in mixtures of platelets 



and ideal polymer 



The 
(DFT) 



undamental-measure theory (FMT) is an (approximate) density-functional theor 
18, originally proposed by Rosenfeld for additive hard sphere mixtures 
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An early extension to convex non-spherical particles has been given in Refs. 
applied to homogeneous and isotropic fluid states, this theory yields the correct second virial 
coefficient of the equation of state, but fails to recover the exact density functional up to 
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second order in density. To remedy the latter problem, an interpolation between the hard 
sphere Rosenfeld functional and the Onsager functional for elongated rods was proposed 



23[ and applied to the bulk isotropic- nematic transition. An FMT that originated from the 



treatment of parel 
model, see e.g. j27 



hard cubes 
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2& was a. 



28l. See also Refs. 



3G, 



so used very successfully for the Zwanzig 



3l| for the treatment of various hard body 



system mixtures. Recently attention was paid to binary mixtures of hard rods and polymers 



32, 
34, 



3c I . Also notable is the treatment of hard body fluids based on two-point measures 

The Bolhuis-Frenkel model of hard spheres and vanishingly thin hard needles [37| can 
be considered as the simplest model hard core mixture of spheres and rods. Previous work 



and an extension to 
Predictions for 



was concerned with the formulation of a DFT for this model 
include rod-rod interaction on the Onsager (second virial) level j3^ 
the orientation orde ring of the rods at a free interface between (isotropic) sphere-rich and 
sphere-poor phases were successfully confirmed by simulations j4l|. Adding a third 
component to this system, the fluid demixing phase behavior of ternary mixtures of spheres, 
rods, and model polymers was found to be rich j^]. Quenching one of the components 
led to investigations of hard spheres immersed in random rod networks 43|, and of the 



isotropic- nematic transition of rods in matrices of immobilized spheres 4J|. The present 
paper is concerned with the introduction of hard platelets to the FMT framework. 

Central to FMT is the so-called deconvolution of the Mayer bond fij between particles of 
species i and j into "weight functions" w'j^ that are characteristic of the shape of particles of 
each species i; the label u indicates the type of weight function. Recall that the Mayer bond is 
related to the pair interaction potential Vij between species i and j via fij = exp(— /JV^j) — 
1, where (3 = l/{kBT), is the Boltzmann constant and T is temperature. For hard 



bodies /j. 



-1 if the two particles overlap and fij = otherwise. In the FMT framework 



a convolution integral is used to express fij as a sum of terms each of which is bilinear 
in the weight functions, i.e. possesses the structure Wi * Wj, where i,j label the species, 
u, T label the type of weight function, and the asterisk indicates the (three-dimensional) 
convolution. The weight functions are characteristic of the geometry of single particles, i.e. 
vanish (already) beyond the physical extent of one particle. Note that this is different from 
the behavior of fij, which extends to larger separations and vanishes provided two particles do 
not overlap. However the convolution of pairs of (shorter ranged) weight functions restores 
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the correct (longer ranged) non-locality of the Mayer bond. This concept was originally 
introduced for hard sphere mixtures 20|. An extension to arbitrarily shaped convex bodies 



was proposed in Refs. [21|, |2j| , which yields the correct second virial coefficients, but gives 
only an approximation for the Mayer bonds. This is insufficient to describe e.g. nematic 
ordering (see e.g. Ref. |22] for a discussion). However, obtaining exact deconvolutions is 
possible in prin ciple, as the cases of the Mayer bond between a sphere and a vanishingly 
thin needle Issl ] and between a pair of thin rods 3^ demonstrate. 

To give an overview of the general structure of FMT, the weight functions wj are used 
to to build so-called weighted densities via convolution with the bare one-body density 
distribution (of each species in the case of a mixture). The weighted densities are then input 
as arguments to an (analytically given) excess free energy density, $. Integrating $ over 
space (and director space in case of anisotropic particles) yields the excess free energy of the 
system, which is the only unknown contribution to the grand potential functional being the 
central quantity in (classical) DFT. For any external potential (like describing gravity or a 
wall) minimizing the grand potential with respect to the one-body density distribution(s) 
gives both the value of the grand potential and the structural information contained in the 
density profile (s). 

Here we show how to treat platelet-shaped particles and mixtures with both spheres and 
needles within the same framework. The relevant weight functions are given and it is shown 
explicitly how all Mayer bonds for the ternary mixture are obtained through convolutions. 
For two-dimensional hard disks the Rosenfeld functional only yields an approximation for 
the exact Mayer bond. Despite this deficieny the two-dimensional Rosenfeld functional 
is considered to be a useful tool to study inhomogeneous situations, see e.g. Ref. Q] for 
an investigation of laser-induced freezing and melting of confined colloidal particles. The 
treatment of the sphere-platelet case inevitably leads to the same defect. We do achieve, 
however, exact deconvolutions of the Mayer bonds for platelet-platelet and platelet-needle 
interactions. 

Furthermore we treat explicitly contributions to the excess free energy functional that 
are of third order in densities. While the treatment of second order terms was a crucial 
input for the original construction of FMT for hard spheres |2(|, this route has been only 
recently pursued in detail 4^, and the authors conclude that they come "close to the edge 
of FMT" in attempting to incorporate properties of the exact expansion of the free energy 
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functional in density on the third virial level into the approximate framework of FMT. In 
particular these delicate cases appear when particle touch, i.e. configurations that are close 
to the edge of the particles. 

Hence we approximate the exact third virial level (the triangle in the diagrammatic Mayer 
/jj-expansion) by terms that are non- vanishing only for cases with common triple intersection 
of the three particles involved. Global prefactors are used to compensate for the "lost cases" 



47| in order to yield reasonable third virial coefficients. 



As a test case for the accuracy of the theory we consider the pure system of hard platelets. 
This is known to undergo a weak first order isotropic-nematic phase transition j3] • In con- 
trast to the application to thin rods, the Onsager (second virial) functional for platelets 
yields only a qualitatively correct account of the phase transition: The transition density, 
density jump, and value of the nematic order parameter in the coexisting nematic are signifi- 
cantly overestimated. Our theory has the same (exact) second order contribution to the free 
energy, but also sports a further third-order term. Higher than third-order terms are absent 



due to the scaled-particle roots pj] of the approach; the vanishing thickness and hence van- 
ishing volume of the platelets effectively truncates the series. Our results show that FMT 
considerably improves on the Onsager treatment, and yields a quantitatively correct picture, 
albeit with slightly too small transitition densities. 

Harnau and Dietrich recently considered bulk and wetting phenomena in a binary mixture 
of colloidal hard spheres and hard platelets an interesting limiting case of the present 
ternary mixture. We postpone a detailed discussion of the relationship of the current work 
to Ref. Q to Sec. HSl 

The paper is organized as follows. In Sec. |H]we define the model of a ternary hard body 
mixture of spheres, needles, and platelets. In Sec. IIIII the Mayer bonds are represented as 
convolutions of (weight) functions characteristic for single particles. In Sec. IIVI we proceed 
to control the third virial level. Readers who are interested in the structure of the resulting 
theory rather than its derivation might wish to skip Sees. IIIII and HVl and directly go to Sec. 
IVl where the weighted densities are defined as convolutions of the weight functions with the 
bare density profiles. Sec. IVII presents the excess free energy functional. Explicit expressions 
for the relevant quantities in planar and uniaxial geometry are given in Sec. I VI II We apply 
the theory to the isotropic-nematic transition of hard platelets in Sec. IVIIII and conclude in 
SecHXl 
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II. MODEL 



We consider a mixture of hard spheres (species S) with diameter a = 2R, where R is the 
radius, hard platelets (species D) of diameter 2i?^ and vanishing thickness, and hard needle- 
like rods (species A^) of length L and vanishing thickness. The pair interaction potential 
Vij between any two particles of species i,j = S,D,N is infinite if their geometrical shapes 
overlap and zero otherwise. The one-body density distribution of species i = S,D,N is 
denoted by P5(r), pi:)(r, fl), and pAf(r, respectively, where r is the position of the particle 
center, and is a unit vector pointing along (normal to) the shape of the needle (platelet) 
describing the particle orientation in space. 

III. DECONVOLUTION OF THE MAYER BONDS 
A. The sphere-sphere and sphere-needle Mayer bonds 

For completeness we first summarize results from the literature for hard spheres and their 
mixtures with needles. Rosenfeld's hard sphere weight functions |2(| are 



where 0(-) is the unit step (Heaviside) function and 5{-) is the Dirac distribution. Here 
and in the following the wl are quantities with the dimension of (length)"^"^; the subscript 
V indicates a vectorial quantity. Further, linearly dependent, weight functions are wf{r) = 
tyf (r)/(47r_R), Wq^t) = tyf (r)/(47r_R^), and wf]^(r) = wf2(i")/(47r-R). For pure hard spheres 
the Mayer bond is obtained through 




6{R-\r\) 
«^f(r)^, 



e(i?-|r|) 



(2) 
(3) 



(1) 



fssir)/2 =wl{r) * w^{r) + w^r) * (r) 
-wf2(r)*w,fi(r). 



(4) 



where the (three-dimensional) convolution is defined as 
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and also implies a scalar product between vectors, as appears in the last term on the right 
hand side of Eq. (jU . 

A pair of vanishingly thin needles does not experience excluded volume interactions, but 
needles do interact with hard spheres. For such a binary mixture the needle weight functions 
used in Ref . ^] were obtained from the prescription of Refs. j2l|, , and are given by 



(5) 



1 /"^/^ 

w;f(r,0) = - / dl5iv-in). (6) 

4 J -L/2 



Introducing a "mixed" weight function that is non- vanishing on the surface of a sphere, but 
carries a dependence on orientation, 

^r(r,0)=2|w,^2(r)-0|, (7) 

allows one to obtain the Mayer bond between sphere and needle through convolutions, 

-fsNir, n) = wliv) * <(r, n) + ^f^(r, O) * <(r, O), (8) 

see Ref. iSftl for the explicit calculation. For the case of residual rod-rod interactions in the 

nn 

Onsager limit, Jnn can also be deconvolved into weight functions, see Ref. laa, m fo^ the 
details. In the following, however, we will restrict ourselves to the case of vanishingly thin 
needles. 



B. Strictly two-dimensional hard disks 

As a prerequisite for our subsequent treatment of platelets in three dimensions, we give 
an overview of Rosenfeld's functional for the model of strictly two-dimensional disks. This 
constitues a prerequisite for considering the overlap between a hard sphere and a hard 
platelet in three dimensions fSec. IIIIEl below). Here we deal with a multicomponent two- 
dimensional mixture of disks with radii Ri of species i, and characterized by the pair potential 
Vij{r) = oo if r < Ri + Rj and zero otherwise, where r is the center-center distance between 
the disks of species i and j. (When regarded from the viewpoint of the three-dimensional 
model, all disks possess the same orienation O, perpendicular to the two-dimensional plane 
of position coordinates; we will not delve into this delicate dimensional crossover in the 
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following). Rosenfeld's weight functions for this model are 



w?{v) = eiR.-\r\), (9) 
w?{T)=6{R,~\r\), (10) 



w2(r)=«;W(r)^, (11) 



r 



and there is an additional linearly dependent weight function, WQ\r) = W2\r)/{2TTRi). The 
exact Mayer bond fij between species i and j is then approximated through fij ^ f*j, with 



(12) 



+ ^ {wf\r) * w'f\r) - w2(r) * w^^V)), 
where here (and only here) * denotes the two-dimensional convolution, 

h{r) * g{v) = J d^x h{x)g(x — v), 

and again the convolution implies a scalar product between vectors. 

As our subsequent treatment of the three-dimensional sphere-platelet overlap relies heav- 
ily on Eq. ((121), calculate the explicit functional dependence of f*j{r) on r and discuss 
some of its properties. We choose two-dimensional Cartesian coordinates r and x, such that 
r = (0, r) and x = {x' simp' ,x' cos ip'), which allows us to write the first term on the r.h.s. 
of Eq. (O as 

= / d<p' x'dx' 6{x' -Ri)Qi R, - \/x''^ + - 2rx' cos ) (13) 

Jo U 27ri?j V ) 



1 .r' + Rf-R^^~ 

— arccos ( 

71 



which is, by the cosine theorem, equal (up to a factor of 27r) to the length of the arc which 
the rim of the disc i traces across the interior of disc j (see Fig.Q). By symmetry, the second 
term on the r.h.s. of Eq. ()12p gives Eq. (fT^ with Ri and Rj interchanged. To calculate the 
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remaining third term on the r.h.s. of Eq. ()12p . we use the same setup as above and obtain 
d'x (w?{^)w['\^ - r) - wS(x) ■ w^^ (x - r)) 

= — y dip'j dx' x'5 {Ri - x') 5 i^Rj - y^r^ + x'"^ - 2rx' cos (/^'^ 



1 

2^ 



^ I ^ ( X smip \ I X simp 

RR' \ J yx' cos ip'—r 



(15) 

^ e{R. - \R, - r\) , r^-iR^~R^r ^__ (ig) 

Armf - (r2 + _ my 



IT 



In the special case that both platelets have the same radius, Ri = Rj, as relevant for 
describing the pure hard disk system, Eq. (fTI)|l further simplifies to 

^ I d'x (^«(x)«;W(x - r) - wS(x) ■ wS(x - r)) = 0(2i?, - r) Z- (17) 

We plot fij{r) in Fig. El for three typical values of the size ratio Ri/Rj = 0.44,0.71,1. 
For small separations r, such that the (small) disk i lies completely inside the (big) disk j, 
the result is exact, /j*(r) = — 1 for r < _Rj — R^. As r increases, —f*j first decreases and 
then increases, exhibiting a divergence at contact, r = Ri + Rj. Beyond this threshold, 
for r > Ri + Rj, f*j is (correctly) indentically zero. Note that the deviations around the 
exact value inside the core, fij = — 1, balance such that the second virial coefficient is exact, 
/rfV/,,(r) = /rfV/*(r) = -7r(i?, + R.f. 



C. Platelet-platelet Mayer bond in three dimensions 

Although the particle shapes appear to be equivalent, overlapping pair configurations are 
very different for freely rotating platelets in three dimensions as compared to the strictly two- 
dimensional case. Note that in three dimensions we only need to consider cases of differing 
orientations of both platelets, 17 ^ fl', as configurations with strictly equal orientations, 
f2 = fl', carry vanishing statistical weight. The intersection volume differs markedly from 
the two-dimensional case (area of intersection), requiring very different weight functions to 
express the Mayer bond. This task can be accomplished exactly, in contrast to the two- 
dimensional case above. The direct application of Rosenfeld's recipe to platelets 
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FIG. 1: Overlapping pair of strictly two-dimensional hard disks with radii Ri and Rj and center- 
center distance r. Shown are the x and y axis of the coordinate system used, the angle a between 
the y-axis and Ri^ as well as half of the part (bold) of the circumference of disk i that lies inside 
disk j. 

yields the scalar weight functions 




(18) 
(19) 




1 



5{R 



r|)5(r-n), 



(20) 
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FIG. 2: Rosenfeld's approximation for the (negative) Mayer bond, —fijir), of two-dimensional hard 
disks of radii Ri and Rj, for size ratios Ri/Rj = 0.44, 0.71, 1 (from left to right) as a function of the 
(scaled) center-to-center distance r/a, where a = 2Rj. In Sec. lIIIEl the same result is obtained for 
the Mayer bond in three dimensions between a sphere of radius R and a platelet of radius as a 
function of the (scaled) distance r/a between the center of the platelet and the axis of the sphere 
(parallel to the orientation fi of the platelet), where a = 2Rj-). The (scaled) distance between the 
center of the sphere and the plane of the platelet is 0,0.7,0.9 (right to left), and R = Rd- The 
result for the Mayer bond obtained is exact in cases (not shown) where the platelet i) fully cuts 
through the sphere or ii) lies completely inside the sphere (provided Ro < R). 

corresponding to the surface, ^2 = 27r_R|,, integral mean curvature, = nR^/A, and Euler 
characteristic, = 1, of the platelets, respectively. The weight function that describes the 
particle volume vanishes, (r, O) = 0, due to the vanishing thickness, and hence vanishing 
volume of the platelets, = 0. 



In contrast to Rosenfeld's treatment 



21 



we here aim at an exact deconvolution of 
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/£)£)(r, f2', n), where r is the center-center distance vector between both particles and ft and 
n' are their orientations. We introduce 



where wf(r, ft), given through Eq. (fTU|) . describes the platelet rim, and ft' is the orientation 
of the second platelet. Keeping Rosenfeld's surface weight function, ti?|'(r, f2) as given in 
Eq. ()18|1 . we recover the Mayer bond between platelets via convolution. 



As the r.h.s. of Eq. consists of two symmetric terms, it is sufficient to consider the first 
one; the second one gives an equivalent contribution to foD- Without loss of generality 
we place the particles such that one platelet is located at the origin with its orientation 
vector pointing up (into the positive 2;-direction), and place the center of the other platelet 
in the y-z plane; see Fig. El for an illustration. The chosen coordinates are: r = {0,r,z), 
n = (0,0,1), ft' = (Vl — sin if, \/\ — cos z\ and the integration variable is x = 
(r' sin r' cos 2;'). Using ()19|) and ()21|) we insert the explicit form ''^(r, fi; f2') = 



^(r, fl; rj') =^|0 • (rj' X r) |wf (r, fl) 



(21) 



/z5D(r, O; U) = w^'^ir, O; W) * (r, U') 

+ w^{r, n)*wf^{r, ft'; ft). 



(22) 
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FIG. 3: Setup for the calculation of the platelet-platelet Mayer bond. Platelet D is located in the 
x-y plane with its center at the origin. Platelet D' has its center in the y-z plane and an arbitrary 
orientation. The calculation checks for the intersection of the rim of D and the interior of D' (bold 
dot). The one-dimensional intersection volume of the two platelets is shown as a bold line. 



13 



\n' ■ (r X f2)| (5 {R^ - |r|) 5 (r ■ fl) /(Si?^) into the first term on the r.h.s. of ^ to obtain 

j d^xw^^{x,n;n')w^{x-v,n') 



2ix 



-'0 J -oo Jo 



dyj' / dz' j dr'r 



W^rr^ain I X ( r' COS ip' 



Vl — 2*^ COS I/P 



(5 (i?^ - + 5 (2') e(i?B - I (r' sin (^', r' cos - r, z' - 2) | ) (23) 



X 6 



r sm (/3 
r' cos ip' — 
z'-z 



\/\~z'^ sin 
VT-S^cos 



27r 



sin (/3 
cos 1/3' 




— \/\—z^ cos (/P 
\J 1— 2*^ sin 




X 0(i?j5 — |(i?^ sin _Rj) COS v?' — r, — z)|) 



(24) 



X 5 



Rj^ sin 
_Rjj cos 93'— r 



-2^ sm 



-2^ cos Lp 



^ [ dip' sm{ip - ip') Q{2rR^cosp>' -r^ - z^) 

2 JO 

X 6 ^rVl — ^2 cos (p — RjjVl — ^2 cos{(p' — (p) + zz 



^ 5Z ® (^ri?^ cos ( ip ± arccos 



RnVT 



^2 



2 2 

r — z 



(25) 



(26) 



X 



6 ( Rj^Vl — -22 — rVl — -22 cos(y9 + 



where the integrations over z' and r' in Eq. (j^Hj) are straightforward. From ()24|) to we 
rewrote the argument of the Dirac delta distribution such that (p' appears only once, using 
the identity sin ip' sin ip + cos ip' cos ip = cos(v9' — ip) . 

Given the complexity of the overlap condition between two arbitrarily oriented platelets in 
three dimensions, it is not surprising that Eq. ()2fij) is an involved expression. It can be viewed 
as counting the number of intersections between the platelet D' and the rim of platelet D 
js^. The second term in Eq. ()22|) counts the number of intersections between platelet D and 
the rim of platelet D'. Hence in total there are two (zero) intersections for an overlapping 
(nonoverlapping) configuration, and indeed the Mayer bond between platelets, Eq. (j221), is 
recovered through convolution of (orientation-dependent) single-particle measures. 



D. The platelet-needle Mayer bond 



We next consider the platelet-needle Mayer bond, /DAf(i", ^d] ^n), where r is the differ- 
ence vector from disk center to needle center, and Qd and rijv are the orientations of the 
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FIG. 4: Setup for the calculation of the platelet-needle Mayer bond. The needle is located at r^r 
in such a way that it points at the origin. The platelet is located at r^i. 

disk and needle, respectively. We define a "mixed" weight function for the platelets that 
is non-vanishing on the platelet surface, but carries an additional dependence on the rod 
orientation, 

^(r, fin; ^n) = 2\ftD ■ n^kf (r, fin), (27) 

where (r, fl), as given in Eq. (fTS|). describes the platelet surface. This allows us to obtain 
the Mayer bond between platelet and needle via 

-/^,v(r, Qd; ^n) = w^'^ir, Qd; ^n) * <(r, Qn). (28) 

The validity of Eq. (j28|) can be seen by employing cylindrical coordinates and putting 
the platelet on the y-axis, its orientation being along the z-axis. Position the nee- 
dle in such a way that its intersection with the x-y plane is the origin (cf. Fig. 
The needle is described by = {zN/zy/l — sin(^, zjsi /z\/l — cos<^, z^^ and fiA? = 
{\/\ — z? sin ip,\/\ — z? cos z). The coordinates of the platelet are r^) = (0,r£),0) and 
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(30) 



= (0,0, 1); the integration variable is x = (a;' sin x' cos 2;'). From ()27|) and ()18p we 
obtain the exphcit expression 102'^ {r, flu; ^In) = 4 6(-R^ — |r|) 6 (r • ri^) ■ ^n\, that 
we insert together with wf^(r, Ojv), as given in Eq. (jH)), into Eq. (PH|) to obtain 

-fDN{r, nn; ^n) = J d'x <(x - r^v, nN)w^^i^ - r^, ^n, n^) (29) 

/oo /•2iT POO /'L/2 

dz' / d(p' / dx'x' / dl 
00 io io J-L/2 

X (5 ^a;' sin — zj^f/zVl — z^ sin + ZVl — sin 
X b ^a;'cosv9' — zn/z^/1 — cos(^ + /Vl — cos(^ 
X (5(z' — 2;^ + /^) Q{Rd — |(x'sin(y9',a;'cosv9' — td, z')\) \z\6{z' 

i*2tt poo 

= / dip' / dx' x'6 {x sin ip') 5 {x' cos cp') 
Jo Jo 

X Q{Rd — |(x'sin(y9',x'cosv9' — r£))\) Q^—\z\ — \zn\^ 

= e{R^-rn)Q(^^\z\-\z^\y (32) 

where the integrals over z' and / are straightforward to calculate. In Eq. (jHT|) . the integral 
over ip' is performed first, rendering the x'-integral trivial. 

The result, Eq. (jHSj), is indeed the Mayer bond between platelet and needle: Due to the 
set-up, an overlap can occur only at the origin and the step functions in Eq. (p?^ provide 
a means to assess whether or not this is the case. The first step function checks whether 
the origin is inside the platelet while the second step function checks whether the needle 
intersects the origin. 



(31) 



E. The sphere-platelet Mayer bond 

Treating the sphere-platelet Mayer bond leads to the most involved geometry of pairwise 
particle overlap in the ternary mixture under consideration. It turns out that this case is 
intimately related to the Mayer bond of the strictly two-dimensional hard disk mixtures. As 
detailed in Sec. lIII Bl the FMT gives only an approximate representation in this case. Still the 

losenfeld functional for two-dimensional hard disks is a reasonably accurate theory, see e.g. 

4^ for a recent study. Here we deal with this problem on the same level of approximation 
as the two-dimensional Rosenfeld case. Note that cutting the sphere with the plane of the 
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platelet yields a disk as the shape of intersection. We need to consider its overlap with a 
platelet, and introduce the weight functions 

w;f^(r, n) = -y/l-{vn/R)^wUr), (33) 

TT 

wff(r,^^) = 4 ^-^^^"^% ^(r), (34) 
wf,(r,^^) = ^<(r,^^), (35) 

where (r), as given in Eq. (0), describes the sphere surface, and (r, f2), as given in Eq. 
(fTTHl . describes the platelet rim. Eq. (fH^ defines a vector field tied to the surface of the sphere 
(with radius R) . In contrast to the "radial hedgehog" of the (classic) vector weight function 
w^2{^)j given in Eq. Q, the direction of the current vector field is radial with respect to 
the (platelet) fl direction and its magnitude decreases towards either pole - a configuration 
one could refer to as a "cylindrical hedgehog" . Note further that wf'^ = |wf^|, in accordance 
with W2 = |w^2l- The vector-field of Eq. (jHSj) resembles a corona and is the straightforward 
generalization of the corresponding two-dimensional hard disk weight function, w^\^(r), given 
in Eq. ((TT)). We use these functions to approximate fsD^ fsD with 

-/|^(r, n) = wiiv) * <(r, n) + wf{v) * (r, Q) 

+ wf^(r, n) * <(r, n) - w,^f (x, n) * w,^i(r, ft), (36) 

where again the convolution implies a scalar product between vectors. The three terms on 
the r.h.s. of Eq. ()36|) represent the arc that the rim of the platelet traces inside the sphere, the 
arc that the sphere traces on the platelet (as illustrated in Figs. El and El respectively), and 
an additional contribution of the cusps where the two arcs meet. We put the platelet into the 
x-y plane with its center at the origin. The center of the sphere is located in the y-z plane, 
see Fig El the coordinates chosen are r = {0,r,z), f2 = (0,0,1), x = {r' simp' , r' cos (p' , z') . 
We obtain for the first term on the r.h.s. of Eq. ()3fi|l 

d^xw^(x,rj)wf(x-r) 



r.27r 



PZTT PCX) POO . 

— — — / d(y9' / d/ / r' d.r' b [Rjj - \fi 
^t^Rd Jo J -00 Jo ^ 



(37) 



X e{R- |x-r|) 

Q{2rR^-y + Rl + z"^ - R^\) . (38) 



1 fr^ + Rl + z^-R^ 

— arccos 

TT V 2rR^ 
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FIG. 5: Setup for the calculation of the sphere-platelet Mayer bond. The platelet is located at the 
origin. The convolution Wq * is equal to the angle under which the arc shown in bold appears. 

By the cosine theorem, this is equal (up to a factor of 2tt) to the length of the arc which the 
rim of the disc traces inside the sphere, see Eq. (jl4j) . Note also that — is the squared 
radius of the circle that the platelet cuts out of the sphere. 



18 




FIG. 6: Setup for the calculation of the sphere-platelet Mayer bond. The platelet is located at the 
origin. The convolution wf * W2 is equal to the angle under which the arc shown in bold appears. 

For the second term on the r.h.s. of Eq. (|36|) we obtain 
j d^xwf(x)wf (x - r, O) 

1 POO p27r poo . 

— / d^' / dy?' / dr'r 5 (r - ^/ r''^ ^ z'A 
A-kR Jo Jo ^ ^ 

X 2 e (^R^ - + r'2 + {z' ~ z)^ ~ 2rr' cos 6 {z' - z) 
^ d<^'Q{R- \z\) e(^RD- -\/r2 + i?2_^2_2rv/;R2^^cos(^'^ 



27r _ 

1 /r2 + R2-z^-R^ 



— arccos , 



D 



e 2rVR^-z^ - y + R^-z^- Rl 



D > 



(39) 

(40) 
(41) 



which is indeed the length of the arc that the sphere traces on the platelet divided by 2tt, 
see Eq (HH). 
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The relevant weight functions to treat the third and the fourth term on the r.h.s. of 
Eq. dnSl) are exphcitly w^^{r,n) = 45(i?-|r|)[r - (r • n)n]/{TiR) and wfi(r,ri) = 
S {R^ - |r|) (5 (r ■ Cl) r/(8i?^), hence 



J d^x [wf^(xX(x - r) - wff (x)wfi(x - r 

■j^ l'2'K roo roo ^ 

= — / dip' / dr' r' / dz'5 ( R — \/ r'^ + z''- 
27r Jo h J -co ^ 

X 6 (^Rf) — A/r^ + — 2rr' cos ip' + {z' — z) 



x6{z' - z) 



r sm(/3 
r' cos tf>' 




r sin \ \ 
r' cos ip'—r 1 



,/2 



RnVr^ 



12 



— / dip' / r' dr'6 [R — V r'^ + z'^ 
X S ( Rd — A/r^ + — 2rr' cos 



„/2 



rr' cos 



2r \/ R? — z'^ cos <y9'^ 



X 



R 

-Q(2rVR 



R 



D 



2_^2 



TT \ 



Rl\] eiR-\z\) 



R^-{VW^ 



^Ar^{R'^ - z^) - (r2 + i?2 _ ^2 _ ^2^)2 



- 0(i?^ - - r|) , ^^^^ 



Rn? 



r^ + R'^- Rlf 



(42) 



(43) 



(44) 



(45) 



(46) 



Considering again that R'^ = R^ — z^ is the squared radius of the circle that the plane of the 
platelet cuts out of the sphere, we recover Rosenfeld's approximation for the Mayer bond of 
hard disks in two dimensions, see Eq. Hence the same defect, as displayed in Fig. |21 

for two-dimensional hard disks, is encountered for the overlap between sphere and platelet. 
Based on the very reasonable quality of the two-dimensional Rosenfeld functional, we expect 
this to do little harm in typical applications of the theory. 
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IV. THIRD VIRIAL LEVEL: TRIANGLE DIAGRAMS 



The contribution of third order in densities to the excess free energy of a mixture has the 
generic exact form of a sum of triangle diagrams, 

-^E/ ^^Z^^^^) / ^^''^^■(^') / dsVfc(s")^(s,s')/,fc(s',s")/.fc(s,s"), (47) 

ijk 

where /c label the species, the sum runs over all species, fij{s,s') is the Mayer bond 
between particles of species i and j with coordinates s and s', respectively, and s (and its 
primed versions) is a composite variable appropriate for the species under consideration, i.e. 
for the case of (uniaxial) rotators like platelets s = (r, fi), J ds = / dr J dri/(47r), and for 
spheres simply s = r, J ds = J dr. 

Even for one-component hard spheres FMT yields only an approximation to (jTTj) due 



to the occurrence of "lost cases" j46l. l47l| : The FMT approximation for the triple product 
of Mayer bonds in (jTfjl is non-vanishing only for configurations r, r', and r", such that 
three spheres with (hard sphere) radius R centered at these positions share a non-vanishing 
common volume of intersection. Note that the existence of such triple intersection implies 
pairwise intersections, but not vice versa: An illustrative example is that constituted by three 



touching spheres; see e.g. 4y| for an in-depth discussion. The lost cases are compensated 
by overcounting the cases with triple intersection, such that a reasonable (exact in the pure 
hard sphere case) third virial coefficient for the bulk fiuid results. The upshot is that new 
weight functions result, that can be, via weighted densities, incorporated to improve the 
accuracy of the density functional approximation. 

Here we follow a similar strategy and aim at controlling the dependence on positions and 
orientations through appropriately constructed weight functions. For triplets of platelets in 
configurations with common triple intersection we can express 

fooir - r',n- Q,')fDD{r - r", Q, ^^")/DD(r - r" , ^^") 

= Addd J dxMif ^^(r - X, n, n'; n")w^{r' - x, n')w^{r" - x, Q"), (48) 

where Addd is a constant, w^(r, 17) is the surface weight function defined in Eq. (fTHj) . and 
we introduce an additional surface weight that carries a dependence on two further (platelet) 
orientations, 

w^^^ir, n, n'; n") = -\n ■ {n' x n")\w^{r, n), (49) 

TT 
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with normalization such that J dr J ^ J ^ J ^^(r, Q'; Q") = = 2ttRI. Set- 

ting Addd = 7r/64 = 0.0491 yields equality in Eq. (PHj) for particle configurations with 
common triple intersection of the three particles shapes. Considering the bulk equation of 
state for pure platelets from SPT, and in particular its (approximative) third virial coeffi- 
cient, the compensation for the lost cases (those without common triple intersection, but 
with pairwise intersections) results in A^ijij = tt/32 = 0.0982. In the following we will 
rather use the value from FMT, Addd = l/(47r) = 0.0796, which is numerically similar 
to the SPT result. The advantage of the FMT scheme is that it consistently applies to 
mixtures, see Ref. js^ for a detailed discussion. 

For the platelet-sphere mixture the triangle diagram stemming from a pair of platelets 
and a single sphere, provided that there is a a common intersection of all three particle 
shapes, can be expressed as 

fsDir; r', n')fnD{r', Q'; r", n")fsD{r; r", Q") 
= j dx[A5Z)^wf^^(r - x; O'; 0")wf (r' - x, 0')w^(r" - x, n") 

+ (r - x) ""(r - X, a- n")w^{r" - x, fJ") 

+ w^{r' - X, n>f ^(r" - X, n", n'))], (50) 



where Asdd is a constant, and we have introduced a sphere weight function that depends 
on two platelet orientations, 

^r^(r; O'; n") = -{{Q' x O") ■ w^r)], (51) 

71 

with the normalization / dr / ^ / ^wf^^(r; 1^") = ^| = AnR\ Again Asdd = 
7r/64 = 0.0491 yields equality in Eq. ()50|) provided that there is common triple intersec- 
tion. We will use the FMT compensation, Asdd = VI^tt) = 0.0796. 

The remaining case is that of triplets consisting of a single platelet and a pair of spheres. 
This case involves similar geometry as the platelet-sphere overlap, and hence we do not aim 
at an exact treatment. Again in the spirit of the two-dimensional hard disk case (cf. Sec. 
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(r' — x)u'f (r" — ^)wq (r — x, f2) 



IIII B|) . we use an approximation, 

fss{r - r")fsD{r'; r, n)fsD{r"; r, ft) 
dx 

+ wlir' - x) (ti;f^(r" - x; f2)^;f (r - x, f2) - wff (r" - x; f2) ■ wfi(r - x, 
+ wlir" - x) («;f^(r' - x; n)w^{r - x, f2) - wff (r' - x; ^^) ■ wf,(r - x, 17)) 
+ AssD w^ir - X, 17) (wi'^iv' - x; 17)«;f^(r" - x, n) 



wff(r"-x;f2).w,^2'^(r"-x;17; 



(52) 



where ^55^1 = 1/16 generates the best possible approximation and Assd = l/(4vr) compen- 
sates for the lost cases in the FMT approximation. 

V. WEIGHTED DENSITIES 



We first summarize the novel weight functions that describe the platelet degrees of free- 
dom in the ternary mixture. From the original Rosenfeld prescription we have obtained [see 
Eqs. (UHll-dlll) in Sec. nTlT!] 



w^{v,n) = 2Q{Rj,-\v\)5{vn), 
<(r,17) = 5(i?^-|r|)5(r-12)/8, 



2txR 



-5(i?^-|r|)<5(r-f^). 



(53) 
(54) 

(55) 



D 



In order to treat platelet-platelet and platelet-needle overlap we have introduced [see Eq. 
da} in Sec. rrrm and Eq. ^ in Sec. mTP] 



w 



f^(r, n') =—\n ■ iSl' X r)|wf (r, 



R 



D 



^(r, fic; 0;v) = 2|ri^ ■ ri;vk^(r, ^d) 



respectively. To treat sphere-platelet overlap we use [see Eq. in Sec. IIII Ej 

4 



^f^(r,fl) 



TT 



1 - (r ■ ^^/i^)2wf (r 



wf2^(r,^2) = 4- ^ ^ 



wfl(r,^^) 



TT 

n 



R 



^f(r), 



R 



<(r,^2), 



(56) 
(57) 

(58) 
(59) 
(60) 



D 
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The consideration of the third order contributions to the excess free energy has led to [see 
Eqs. (gni) and dSU in Sec. |IV| 



f ^^(r, n, n'; n") = -\n ■ {n' x n")\w^{r, n), (ei) 



7r 

wf^^(r; n'; Q") = x Q") ■ wf,(r)|. (62) 

TT 

All weight functions are used to build weighted densities via convolutions with the bare 
density profiles. For hard spheres the scalar and vectorial weighted densities are respectively 
given by 

nUr) = wUr)*Psir), ^. = 0,1,2,3, (63) 
n5(r) = w^(r)*p5(r), u = vl,v2, (64) 

where ^^(r) and w5(r) as given in and below Eqs. CQ)-© are the classic Rosenfeld hard 
sphere functions. Adding non-spherical particles requires convolving the sphere density 
profile with orientation-dependent weight functions. The coupling to needles is accomplished 
via 

nf^(r,f7)=«;f^(r,f7)*p^(r), (65) 

where wf^(r, Q) as given in Eq. ((7j) is a mixed sphere-needle function. Treating the platelets 
follows a similar outline, 

nf^(r,n)=ti;f^(r,f2)*p5(r), (66) 
n,^f(r,f7)=wf,^(r,f2)*p5(r), (67) 
nf^^(r; n, ^2') = ^f^^(r; n, ^2') * ps{r), (68) 

with W2^{r, ft) and ^^^{r, fl) given in Eqs. and respectively, and wf^^(r; ft'; fl") 
as given in Eq. (j^D). 

For needles the pure weight functions are 

(r, n) = <(r, n) * pjv(r, fl), z/ = 0, 1, (69) 

where WQ(r,Q) and w^{r,Q) are given in Eqs. (0) and (jH)), respectively. In case one 
considers needle-needle interactions, additional weighted densities need to be taken into 

rn 

account [39.] . 
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The degrees of freedom of the platelets are modelled through the weighted densities 

n^(r, n) = wj^ir, n) * poir, fl), z/ = 0, 1, 2, (70) 
nfi(r, n) = wfi(r, O) * pz)(r, n), (71) 

^(r, n') = I (r, O; n') * p^(r, ^^), (72) 

^^(r; fl') = J ^^^^(r, fl"; f^; 17') * p^(r, O"), (73) 

with the weight functions w^(r, O), (r, fi), and Wq^t, ft), as given in Eqs. (fTHj). (dH), and 
(120)), respectively, and w^^{r,ft) as given in Eq. The double orientation-dependent 

weight function ^(r, fl; f2') is given in Eq. (|2T| . The weight function that we have 
introduced to control the triangle diagrams, W2^^{r, O'; ft"; O'"), is given in Eq. ()49j) . 
The remaining coupling of the platelets to the needles is achieved via 

nr(r, OO = I ^^r(r, n, n') * p^(r, O), (74) 

where we have used W2^{r, fl) as given in Eq. ()27|) to build a mixed weighted density. 



VI. FREE ENERGY FUNCTIONAL 



In Ref. 38| the excess Helmholtz free energy functional for particles with rotational 



degrees of freedom was expressed not only as a spatial integral (as was proposed in Refs. 
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221), but also as an integral over the director space. In our present explicit consideration 



of three particle correlations an additional integration over director space is required, 

= /''"■/ ^ / 

where the (reduced) free energy density, $, is a function of the set of weighted densities, 
{nl}, where i labels the species and z/ labels the type of weight function. 
For hard spheres $ = $s with 



<l'5 = -n^ln(l-nf) + 



1 — 



|(nf )3 - nfnga • 
^ 87r(l-nf)2 ' ^""^ 

where we have suppressed here and in the following the dependence of the weighted densities 
on the space coordinate r. Eq. (f7B|) is the original Rosenfeld form {2^ that yields, for bulk 
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fluids, the (reduced) free energy density as obtained by tlie Percus-Yevick compressibility 
or scaled-particle) route. Improved versions exist that feature exact dimensional crossover 



51 



52l ]. and the Carnahan-Starling equation for bulk fluid states js^. 
For binary mixtures of spheres and needles the reduced free energy density is $ = $5 + 
$2JV, where $5 is given in (j76p . and ^sn describes the effect of needle-sphere interactions 
38[, given as 

^SN = -n^m ln(l - nl) + (77) 

i rig 

where we suppress the dependence of the weighted densities on r, but make the dependence 
on orientation explicit in the notation. The free energy for bulk fluid states obtained from 



dZZj) is 

theory 



;he same as that obtained from a thermodynamic (free-volume-like) perturbation 
that results in a fluid-fluid demixing binodal that compares well with results 



from simulations. 

For a system of pure platelets $ = with 

$^ = nf ^(^^)n?(0) + (f^)n?^^(f2, n')n^{n'), (78) 

where the prefactor of the third order term in density is according to Rosenfeld A = 
l/(247r) = 0.01326. Scaled-particle theory 0, Q gives a slightly different value, namely 
A = 7r/192 = 0.01636; see Ref. js^ for a detailed comparison of the predictions of these 
approaches for bulk fluid properties. Eq. (fTSj) is the simplest form of a free energy density 
that features the correct second virial level and a similar free energy for bulk isotropic fluid 
states as scaled-particle theory. 

For binary sphere-platelet mixtures we obtain $ = $5 + where the hard sphere 

contribution $5- is given in Eq. (fTHjl . and 

^SD = -n^in)\nil-nl) 

^ nfnf(0)+nf^(l^)nf(0)-ngf(0)-nf,(») ^ n?^{n)n^{n) 

1 — raf 1 — ^3 

nl^mnl^in) - nff (17) ■ nff (17))nf (17) 

+ 



87r(l - nl 



S\2 



87r(l-r2f)2 ^ 247r(l-r2f)2 " 

For a binary platelet-needle mixture the free energy density is $ = + ^dn-, where 
the pure platelet contribution $£) is given in Eq. (|75|) . and a simple second-order coupling 
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between both species is introduced via 



$^^ = n^^(^^;^i>f(^2'). (so) 

Due to the vanishing thickness of particles of both species, any higher order terms vanish 
according to SPT. 

Finally, for ternary sphere-platelet-needle mixtures the free energy density is 

where the first three terms on the r.h.s. are given through Eqs. (|76|) . (|77jl and (|79|l . respec- 
tively. 

VII. PLANAR GEOMETRY 

In many practical situations one is faced with inhomogeneities that depend only on a 
single (Cartesian) coordinate, while being translationally invariant in the two remaining 
directions. A smooth planar wall, where the fluid density profile(s) only depend on the 
perpendicular distance, from the wall is a primary example. For rotators an additional 
simplification arises from (cylindrical) rotational symmetry around the z-axis. Hence such 
problems are characterized solely by z, and the tilt angle d between the particle orientation 
and the z-axis. We give in the following explicit expressions for (reduced) weight functions 
appropriate for efficient numerical treatment of such situations. For completeness we also 
give the results for spheres and needles Soj ]. 

Let us start by defining the full set of reduced densities; we use the parametrizations 
r = (x, y, z) and = (-(9, </)). For the spheres 

wl{z) = JdxJ dyw^^iv), i^ = 3,2,l,0,v2,vl. (82) 

Further mixed weight functions couple spheres to needles, 

wl^{z,^) = Jdxj dy^f^(r,0), (83) 
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and spheres to platelets 

wl^'iz^d) = JdxJ dywr{r,n), (84) 

w.T(^,^) = /^/d^/ dywff (r,f2), (85) 

^"'^'^'^ ^/^/^/^^/ dl/^f^^(r,^^,^^'), (86) 



For the needles 



w^{z,^) = JdxJ dyw^{r,n), u^O,l, (87) 

For the platelets the effective weight functions for planar geometry are obtained by inte- 
grating over the lateral coordinates 

w^{z,i9)^ JdxJ dyw^{r,n), i^- 0,1,2, (88) 



and 



<iz,^) =1^1^^ I dy<ir,n), (89) 

w^''{z,^,&)^ J ^ Jdx J dt/<^(r,0,f^') (90) 

w^^^iz,^,^',^") =J^J^JdxJ dyw?^^{r,n,n',n") (91) 

w^''{z,^n;^N) ^ J^JdxJ dyw^''{r,nn;n^). (92) 
Using the effective weight functions the weighted densities can then be written as 

n'^,^iz,^)^w'^,^iz,^)*ps{z), (93) 



n^{z,^)=w^{z,{})*pry{z,^), u = 0,1,2, (94) 

n^,{z,i^)=w^,{z,^)*PD{z,i^), (95) 

nf^{z,'d')^l [ d^s\u^w^^{z,'&;'&')* pd{z,'&), (96) 
2 Jq 

{z,'&,d')^\ r d^" {z,r;'&;d')*pD{z, ^ ) (97) 

nf^{z,^') = - [ d^ smdw^^ {z,&,d') * pD{z,d), (98) 
2 Jo 
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where in Eqs. fl93p - ()98p (and only here) * denotes the one- dimensional convolution, g{z) * 
h{z) = J g{z' — z)h{z'). In contrast to the needle case the integrands depend non-trivially 
on (f. 

We finally give all results for the effective weight functions in planar geometry. For hard 
spheres the well-known results are 



wliz) 
wl{z) 



-niR^ - z^)Q{R- \z\), 

2'kRQ{R-\z\), 

2^zQ{R-\z\)e^, 



(99) 
(100) 
(101) 



where is the unit vector pointing along the 2;-axis. The linearly dependent weight functions 
are wf{z) = 0(i? - \z\)/2, w^{z) = e{R - \z\)/{2R), w^^{z) = zQ{R - \z\)eJ{2R). The 
mixed sphere-needle weight function is given through 



wi^{z,^) = { 



8Vi?2 sin2^9-z2 

' 2COt(l?) 



X arcsm 

4:'k\z\ cosd 


Further mixed sphere weight functions are 



if \z\ < Rsini) 

if /?sin?9 < 1^1 < i? 

otherwise. 



(102) 



M;f^(2,^9) = Jle(i?- 1^1) / difijR^- (R^ - z^)cosipsm^ + ZCOS^ 



2n 



w,T(^,^9) 



57r 



R 



z sin^ Rd - z^e{R- \z\)e. 



D 



W. 



SDD 



{z,d,d') = Q{R-\z\) / d^' 
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if e > 



i| ^ — f^ + t arcsin ^ j else, 



(103) 
(104) 

(105) 



where is a unit vector in the ^-direction, t = 2 sin sin -(9' sin y?', and s 
^R? — z'^\J sin^ ■{} cos^ ^' + cos^ "i? sin^ ^' — 2 sin ^cos^ sin -d' cos "d' cos . 
For the needles the effective weight functions are 

L 



w^{z,d) = (4cost9)"^e(^-cos^9- l^l); 
w^{z,^) 



-§( — cos?? — \z\ 
2 \2 ' 



(106) 
(107) 
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Carrying out the integrations in Eq. ()88p yields 



' Vi?^sin^92_^2 ^ ^ 

hVRd smw"' — 

Comparing Eqs. p09|) and piO|) to the corresponding expressions for needles, Eqs. ()10fj|l 
and ()l()7j] . reveals that the different particle geometries imprint marked differences in the 
functional dependence on z and Clearly this can be traced back to the differences in the 
full weight functions for needles [Eqs. (0) and (0)] and platelets [Eqs. flT^ and (j^ ]. Note 
that for the present case of vanishing particle thicknesses, w^(^, has no needle counterpart 
to compare with, as the needle surface vanishes. 



Rr,^ sin^ — z'^ 



wfi(^,^) = ^..a I V; e., (Ill) 



8Rj^ sin^ ^9 V _ ^2 



w^''{z,^,r) 



if p2 < 



^ cos t?' 



if p2 sin(^9 - ^9') sin(?9 + > z^ sin^ d' 



y2 \ m2\ sin^ 



, ,, Z'^ + V 
7rp V V ^ ' sin - 



+p| cost?' I 



X arcsm ■ 



p| cot ■d'\ 



else, 



(112) 



where p = a/ -R|, sin^ d — z"^. 



w^^^{z,d,d\d") = -JrI--^. r d^"u{^,^',^",^")e{Rlsin'^ - z') (113) 
TT^ V sm V In 
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where 



cos 'd' sin sin 
2 ' >» sin2 



if > 1 



sin cos ?9 cos y?" — sin cos t?" 

1/2 



7r| sin i9| 

+ sin('^' + ^) sin(t9' - ^) sin^ t?" sin^ 



+ sin^^sin 1?" cos??' sin (y?" arcsin -(9', (yj") j else, 



;ii4) 



and 



(sin t?" cos COS ^ — sin cos i?")^ + cos^ sin^ {}" sin^ 

( 

For the remaining mixed weight functions we obtain 



sin -t} sin -d" cos ^' sin ip" 



sin 



115) 



4 cot {}d cos'dNvR 



e{RjjSin§D-\z\) . (116) 



VIII. RESULTS 



Aheady the binary subsystems with platelets (leaving aside the full ternary mixture), 
i.e. the sphere-platelet and rod-platelet mixtures, are expected to display rich bulk and 
interfacial properties due to the competition between depletion and orientational order. In 
order to find a simple yet demanding test case for the present theory, we restrict ourselves 
further to the pure system of platelets. In bulk this undergoes an isotropic-nematic phase 
transition, where poir^Q,) = po = const in the isotropic phase, and p£)(r, il) = Pd{G) in 
the nematic phase being peaked around the nematic director that indicates the preferred 
direction of alignment of the particles; 9 is the angle between VL and the nematic director. 
The strength of the nematic order is conveniently measured via the nematic order parameter, 

S = Pd j ^PDmP^icose), (117) 

where P2{x) = (3x^ — l)/2 is the second Legendre polynomial. 

We use the generic form of the free energy, Eq. (f75j) . as a multiple integral over a free 
energy density appropriate for pure platelets, $£> as given in Eq. (f7S|). This function depends 
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on the weighted densities nf*, nf*^, and n^-^-^, as defined in Eqs. (f7n|) (with z/ = 2), (f?^ . and 
(f7S|) . respectively. The weighted densities in turn are buih as convolutions of the bare platelet 
one-body density, p(r, fl), with the weight functions W2, wf^ , and W2^^, as given in (|18|). 
(|2ip. and (j49j) [and summarized in Eqs. (j53j) . (j56j) . and respectively. To carry out the 

actual calculations, we have started from the corresponding expressions in planar geometry 
(see Sec. IVII|) . that we have further (numerically) reduced to ^-independent quantities, as 
appropriate for the spatially homogeneous isotropic and nematic phases. 

Alternatively, our FMT excess free energy functional for hard platelets can be viewed as 
being composed of the exact second-order virial contribution and a further (approximate) 
third order term that i) vanishes for configurations without common triple intersection of 
three platelets, and ii) is constant for configurations with common triple intersection; the 
constant is adjusted to give a reasonable value for the third virial coefficient. Note that 
the theory for the full ternary mixture has a similar structure, but features higher order (in 
density) contributions that arise from the finite volume, and hence finite packing fraction, 
of the hard spheres. 

We minimize the grand potential functional 



where /id is the chemical potential, the dependence on temperature and volume has been 
suppressed in the notation, and we formally set the thermal wavelength equal to the platelet 
radius, A/j = R^, to fix an arbitrary additive constant to the chemical potential. For a 
given value of the (scaled) chemical potential fi* = (3 fin, we minimized fl with respect to 
P£)(f2) in either phase. Inserting the resulting distribution into Eq. piSj) yields the grand 
potential of the system, and the condition for equahty of the grand potential in both phases 
locates the phase transition. 

The numerical implementation uses free minimization, i.e. no parametrized form of p_d(^) 
is assumed a priori. We use an equidistant grid in ^-space with 100 grid points in the interval 
[0, vr/2]. The integration over y?" in Eq. (jll3|) to obtain W2^^ was done with 200 grid points 
in [0, 2n]] this only needs to be performed once at the start of the calculation. We have used 




(118) 



where the ideal gas contribution is given by 




(119) 
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FMT 


simulation 


Onsager 


„iso p3 


0.418 


0.473 


0.667 


„nem p3 

Pd -"-d 


0.460 


0.509 


0.845 




0.041 


0.036 


0.178 




0.492 


0.37 


0.781 



TABLE I: Results for the bulk isotropic-nematic phase transition from different approaches. Given 
are values for the (scaled) coexisting densities in the isotropic and nematic phase, p^^R^ and 
P^™R\), respectively, the (scaled) density jump, ApjjR^ = (pj)™ — Pd°)-^D' value of the 

nematic order parameter in the coexisting nematic phase, Snem, as obtained from FMT, and the 
simulation results by Frenkel and Eppenga [4^, as well as theoretical results obtained from the 
Onsager second virial functional for hard platelets. 

alternatively fixed-point iteration or simulated annealing to minimize the density profile; 
results from both approaches were consistent. Up to 1000 steps were sufficient to obtain 
convergence. 

The resulting values of the densities in the coexisting isotropic and nematic phases, 
and p^™, respectively, as well as the value of the order parameter in the coexisting nematic 
)hase, S'nem, are in good agreement with Frenkel and Eppenga's classic simulation results 
3|, see Tab. QJ The FMT values somewhat underestimate the coexistence densities and 
give a too high value for S'nem- The Onsager second virial functional, on the other hand, 
overestimates significantly the coexistence densities as well as S'nem- Hence the actual phase 
transition, as characterized by the jump in coexistence densities and in order parameter, is 
weaker than predicted by the Onsager treatment. To further illustrate our findings, we plot 
in Fig. [7|the variation of S with fi* and in Fig. |Hlthe variation of S with PdRd- We also plot 
the behavior of the metastable nematic branch, i.e. for statepoints where the nematic phase 
remains locally stable against small fluctuations, although globally the grand potential for 
the isotropic phase is the lower one. It is apparent that the value of S'nem very sensitively 
depends on the precise location of the phase transition. Taking this into account we find 
the overall the agreement of the FMT results with simulation data to be very reasonable, 
and to give confidence both into applications to interfacial phenomena of platelets, as well 
as into the accuracy of the FMT for the full ternary mixture. 
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FIG. 7: Nematic order parameter S for a bulk system of pure hard platelets as a function of 
the (scaled) chemical potential, /i*. Shown are results from FMT in the stable (solid line) and 
metastable (dashed line) nematic. Values in the coexisting nematic and isotropic phases, /Snem and 
'S'iso = 0) are shown as obtained from FMT (squares) and the Onsager-type second virial treatment 
(dots). 

IX. CONCLUSIONS 



In conclusion we have derived a geometry-based density functional theory for hard body 
mixtures of spheres, platelets and needles. Both the needles and the (circular) platelets pos- 
sess vanishing thickness and hence constitute the simplest examples of prolate and oblate 
model particles, respectively. Our treatment of the mixture is based on the so-called decon- 
volution of the Mayer bond into single-particle functions which vanish beyond the extent of 
the particle. The Mayer bonds are recovered upon convolution of the single-particle func- 
tions. The construction of the full functional relies further on controlling the third virial 
level and on Rosenfeld's scaled-particle treatment. In order to facilitate future applications. 
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FIG. 8: Same as Fig. [3 but as a function of the (scaled) density of platelets pdR^,- Note the 
considerable reduction of the coexistence density gap obtained from FMT as compared to the 
Onsager-like treatment. 

like wetting of planar walls or capillary phenomena in planar slits, we have given explicit 
reduced expressions for the relevant quantities in planar geometry. 

In a recent contribution Harnau and Dietrich propose and apply a DFT for binary platelet- 
3|. The y o btain a platelet-sphere functional by starting from the rod- 



sphere mixtures 



sphere functional of Refs. 
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and in particular from the explicit expressions for the 
needle weight functions in planar and uniaxial geometry Modifying the definition of the 
relevant angle between the particle orientation and the z-axis the rod weight functions 
are taken to play the role of platelet weight functions. The resulting excess free energy 
functional is linear in the platelet density, limiting the theory to small densities of platelets. 
Recall that for binary mixtures where one component (the depletant) is ideal, the absence 
of higher than linear order terms in the density distribution of this component is a good 
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approximation. Examples are the above rod-sphere mixture and the Asakura-Oosawa model 
of colloid-polymer mixtures, where the polymers are described as non-interacting spheres. 
While in these cases the pure depletant system is an ideal gas, pure platelets constitute 
an interacting system. Hence higher than linear order terms in the platelet density should 
only be irrelevant at low platelet densities. The current theory is in accordance with scaled- 
particle theory 5^ and yields second and third-order contributions to the excess free energy, 
but higher order terms are absent. This behavior can be traced back to the vanishing volume 
of the platelets. 

In the present work we have shown explicitly how our DFT for the sphere-platelet mixture 



• 3, 



we have 



reduces to the correct low-density limit. In contrast to the procedure in Ref. 
treated the full three-dimensional problem and have obtained, up to a defect already present 
for two-dimensional hard disks, the deconvolution of the sphere-platelet Mayer bond and 
hence the appropriate platelet weight functions. Projecting those to planar and uniaxial 
symmetry (appropriate for fluid states at a planar smooth wall like investigated in Ref. 
reveals that the expresssions differ markedly from those for sphere-rod mixtures. This might 
come as no surprise given the fact that the genuine shapes of the particles are one of the 
building blocks of the geometry-based DFT. An immediate consequence is that the platelet 
weight functions, and hence the form of the excess free energy functional, also differ from 
the expressions used in Ref. 

Our successful description of the bulk isotropic-nematic transition of hard platelets can 
be viewed as the first genuine FMT treatment of liquid crystalline ordering in a contin- 
uum model (without any interpolation as e.g. in Ref. [2^). Possible future applications of 
our theory include capillary and wetting phenomena, influence of gravity or other external 
fields, and the study of free interfaces between demixed (and possibly liquid crystalline) 
phases. Further testing the accuracy of the theory, i.e. for the predictions of the bulk di- 
rect correlation functions, see e.g. [^, against results from computer simulations is clearly 
desirable. 
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